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ABSTRACT 

Context. Dust emission is increasingly used as a tracer of the mass in the interstellar medium. With the combination of Planck 
and Herschel observatories, we now have both the spectral coverage and the angular resolution required to observe dense and cold 
molecular clouds. However, as these clouds are optically thick at short wavelengths but optically thin at long wavelengths, it is tricky 
to conclude anything about dust properties without a proper treatment of the radiative transfer 

Aims. Our aim is to disentangle the effects of radiative transfer and dust properties on the variations in the dust emission at long 
wavelengths. This enables us to provide observers with tools to analyse the dust emission arising from dense clouds. 
Methods. We model cylindrical clouds with visual extinctions between 1 and 20 magnitudes, illuminated by the standard interstellar 
radiation field, and carry out full radiative transfer calculations using a Monte Carlo code. Dust temperatures are solved using the 
DustEM code for amorphous carbons and silicates representative of dust at high Galactic latitude (DHGL), carbon and silicate grains 
coated with carbon mantles, and mixed aggregates of carbon and silicate. We also allow for variations in the optical properties of 
the grains with wavelength and temperature. We determine observed colour temperatures, r,.o;o,„-, and emissivity spectral indices, 
ficoioiir, by fitting the dust emission with modified blackbodies using a standard fitting method, in order to compare our models with 
observational results. 

Results. Radiative transfer effects can explain neither the low Tcohm, the increased submillimetre emissivity measured at the centre 
of dense clouds, nor the observed Pcoiour - Tcoiom- anti-correlation for the models considered. Adding realistic noise to the modelled 
data, we show that it is unlikely to be the only explanation of the pcoiour - Tcoiour anti-correlation observed in starless clouds, which 
may instead be explained by intrinsic variations in the grain optical properties with temperature. Similarly the higher submillimetre 
emissivity and the low Tcoiom- have to originate in variations in the grain optical properties, probably caused by their growth to form 
porous aggregates. We find it difficult to infer the nature of the grains from the spectral variations in their emission, owing to radiative 
transfer effects for A < 300 /im, and to the mixture of different grain populations for longer wavelengths. Finally, the column density 
is underestimated when determined with modified blackbody fitting because of the discrepancy between Tcoiour and the "true" dust 
temperature in the innermost layers of the clouds. 

Key words. Radiative transfer - ISM: general - ISM: clouds - Dust, extinction - Submillimeter: ISM - Infrared: ISM 



1. Introduction 

Even if the main lines are understood, a number of questions re- 
main regarding the processes that form a star from a cloud of gas 
and dust. Recent observations indicate that the Galactic interstel- 
lar medium (ISM) is organised as a network of elongated fila- 



clouds, through the emission of the cold dust they contain. Much 
can be learned about the cold clumps from this emission as 
the properties of the grains are known to vary strongly from 
the diffuse ISM to the centre of dense molecular clouds (Tab. 
1). For instance, the temperature of the grains decreases from 
the diffuse to the dense medium (|Schnee et al. 20101, while 



i"'"? Hi^r"';^.''"'''"'' n ^M2010| [Andre et al.|2010[g^ their submillim etre e missivity incre ases (|Cambr6sy et aL 



|et al.|2010 1 and that pre-stellar cores seem to form preferentially 
inside the densest of these filaments (Hill et al. 201JJ |Nguyen| 
|Luong et al.|2011) l. However, the details of either the formation 
of the filaments or their later fragmentation into dense molecular 
clumps are still poorly known. 

The physical, chemical, and dynamical state of the filaments 
is expected to influence the star formation efficiency and the 
mass of the new-born stars. It is thus crucial to characterize 
the properties of the gas and in particular, in the context of 
the Planck and Herschel missions, of the dust grains that com- 
pose Galactic interstellar filaments. Indeed the combination of 
these two observatories offers both the spectral coverage and 
the angular resolution necessary to enable the detection and the 
characterisation of dense, molecular, and potentially pre-stellar 



Stepnik et al.|[2003| IRidderstad et al.||2006| [Kiss et al. 



2001 



2006 



Lehtinen et ai.ll2007! 'Schn ee et al.||2608t [Juvela eraL]|20TT 
Planck Collaboration 2011d]l. There is also observational evi 
dence that the small grains disappear towards dense molecular 
clouds, as the deficit of emission in the mid-IR ( [Bernard et al.| 
1999; Stepnik et al. 2003; Flage y et al.||20091 l. In addition, the 
ratio of total to selective extinction Ry increases from a stan- 
dard value of 3.1 in the difiiise medium to up to ~ 5.5 towards 
dense clouds ([Fitzpatrick & Ma ssa 1988 i [Mathis & Whiffen] 
T989t|Cardelli & Clayton] l"99TT]Campeggio et al.|2007| l, which 



can be explained by the disappearance of the small grains ( Kim 
|et al.|199"4l ). This indicates that the grains grow from the diffuse 
medium to the centre of dense molecular clouds. It also appears 



1 



N. Ysard et al.: Modelling the dust emission from dense interstellar clouds 



to be increasingly clear, though still debated, that dust emission 
at long wavelengths is more complex than the usual modified 
blackbody at temperature T, with a constant emissivity spectral 
index /3 ~ 2 (Tab. 1). Observations towards the diffuse gas as 
well as the dense medium show an anti-correlation between /3 
and T ([Pupae et al.|2003[ [Desert et al.,200 8 ; P aradis etal.|2010 



[Veneziani et al.[[20i0[ [Planck CollaboiiHon||2011d| c'). This be 
haviour can be explained by models of the physics of the emis- 
sion of amorphous solids ( Meny et al.!2007 ; Paradis et al. 201 l| l, 
and is also corroborated by laboratory measurements on inter- 



stellar dust analogues of amorphous carbons (Mennella et al. 



1998 ) and amorphous silicates ( Agladze et al. 1996 ; Boudet et al. 



2005,,Coupeaud et al. 2011 j . ,Malinen et al.. ( .20ir ) also showed 



that a/3 - T anti-correlation can be produced by a warm source 
inside the clouds. In addition, Shetty et al. (20091 found that 
this can naturally result from noise uncertainties and suggested 
that the anti-correlation may also be produced by temperature 
variations along the line-of-sight. Furthermore, laboratory mea- 
surements and models also predict spectral variations in the dust 
emissivity ([Boudet et al.|2005[|Meny et al.|2007[ [Mennella et al. 



1998[[Coupeaud et al. 201 l[)7which have already been observed 
m the Galactic ISM ( [Paradis et al.[[2009) [PlaiTck Collabora"tion[ 
[201 Idp . 

Thanks to the new generation of instruments available, and 
in particular to the combination of Planck and Herschel data, we 
are now able to observe dense and cold molecular clouds, and to 
test detailed models of grain emission and growth. However, to 
describe these dense media, where the illuminated edge is much 
warmer than the shielded centre, detailed calculations of the 
radiative transfer are compulsory. The dust emission spectrum 
arising from dense regions has already been studied by various 
authors. For instance, Fischera & Dopita ( 2008[ l and Fischera 



only radiative transfer and/or noise in the data? What intrinsic 
relation would be needed to reproduce the observed 13 - T anti- 
correlation, when the effects of radiative transfer are considered? 
What is the impact of the variations in the grain properties and of 
the radiative transfer on the quantities that are usually measured: 
column density and dust opacity? These are the main questions 
that we try to answer in this paper. Our aim is to give key in- 
formation to the observers to be able to analyse quantitatively 
observations of dense starless clouds. 

The paper is organised as follows. In Section[2] we describe 
the dust emission model and the radiative transfer code we use 
throughout the paper. In Section [3] we detail the characteristics 
of the dust emission from dense interstellar filaments when as- 
suming that the grains have the same optical properties as in the 
high Galactic latitude diffuse ISM ( Compiegne et al.|201 1 1. We 
explore the effects of both radiative transfer and noise on the re- 
lation between the emissivity spectral index and the colour tem- 
perature. In Section |4] we present the effects of grain growth 
on the emerging dust emission. We consider two scenarios of 
growth, accretion and aggregation, and investigate the relation 
between the colour temperature and the grain physical tempera- 
ture. In Section [5] we show the impact of intrinsic variations in 
the dust optical properties with wavelength and then with tem- 
perature. In Section [6] we test the accuracy of the estimate of 
the hydrogen column density by fitting the dust emission with 
a modified blackbody and investigate how the variations in the 
dust emissivity from diffuse to dense regions can be used to 
probe grain growth. Finally, we present in Section [7] our con- 
clusions. 



2. Models and methods 



( [201 1) investigated the link between the dust emission spectrum 2.1. Dust emission model 



in spherical, passively heated cores and their environment (ra 
diation field, external pressure). They showed that the radiative 
transfer effects lead to an overestimate of the sizes of the cores 
and to an underestimate of their masses, in agre ement with[Evans 



Etar('2Q0r);'Stam atellos & Whitworth[ ( [2055] l; [Stamatellos et al. 



(2004) and Malinen et al. (2011 1. These authors used various 



dust models and cloud geometries, either spheres or cylinders, 
and either symmetric or not, and considered them as embed- 
ded or directly submitted to the interstellar radiation field, but 
all concluded that the radiative transfer has a very strong im- 
pact on the observed dust emission. They showed that the shape 
of the emission profiles differs from the shape of the column 
density profiles, that the apparent temperature diverges from the 
dust physical temperature, and also that the sizes and the masses 
of the clouds cannot be accurately estimated by fitting the dust 
emission with a modified blackbody (Stamatellos et al.,,2010, 
Fischera 2011 Malinen et al. 201 1\ . All these results empha- 



sized the importance of a proper treatment of the radiative trans- 
fer if one wants to trace the variations in the dust properties from 
the diffuse to the dense medium, or otherwise the importance of 
using suitable corrections. 

Therefore, the first question to answer concerns our ability 
to disentangle the effects of radiative transfer, i.e. the mixing of 
different dust temperatures along one line-of-sight, and the ef- 
fects of the variations in the grain properties on the dust emission 
coming from dense clouds. Can we link the "true" dust temper- 
ature at the centre of the clouds to the average temperature we 
measure? How should the grains grow to have the low temper- 
atures and the higher submillimetre emissivity they exhibit in 
dense clouds? Does the observed /3 - T anti-correlation origi- 
nate from intrinsic properties of the grains or is it an effect of 



We use the dust model, DustEIVf] described in'Compiegn e et al 
( |20TT) . DustEM is a versatile numerical tool that computes dust 
emission and extinction as a function of the grain size distribu- 
tion and optical properties. Our standard dust model is defined as 
the DustEM model for dust in the diffuse ISM at high Galactic 
latitude (DHGL according to the nomenclature of Compiegne[ 
et al. ( 201 l| l). Three dust populations are taken into account: 
interstellar polycyclic aromatic hydrocarbons (hereafter PAHs), 
amorphous carbons (small, hereafter SamC, and large, LamC), 
and the so-called astronomical silicates (hereafter aSil). The op- 
tical properties of the SamC and LamC grains come from the 
laboratory measurements of jColangeli et al. (19951, who mea- 
sured the opacity of a sample of sp^-nch hydrogenated amor- 
phous carbon for A < 1900 fim ( Compiegne et al.][2011 i For 
longer wavelengths, the absorption efficiency is extrapolated us- 
ing a single emissivity spectral index j3 - 1.55, equal to its av- 
erage value for 800 < /I < 1900 jum. The optical properties of 
the aSil population were chosen to reproduce the observations 
of the diffuse ISM ( Draine & Lee|1984 ) and have a single emis- 
sivity spectral index equal to 2.11 in the far-IR and in the sub- 
millimetr^ Then, for the small grains (PAHs & SamC), we use 
log-normal size distributions with central radii aq and widths cr. 
For the larger grains, we use a power-law distribution a", starting 
at a = 4 nm, with an exponential cut-off e^^^"^"'^^"'^^ for a > a, 
(Weingartner & Draine 2001 ; Compie gne et al.|201 1 1. The dust 



' Available at htpp : //www . ias . u-psud . fr/DUSTEM 

- We do not use the update of the optical properties made by 



Draine 



( 2003b I, who added a flattening of the emissivity spectral index in the 
submillimetre (p = 2.\ \ — > 1.7). Spectral variations in p are addressed 



separately in Section 5.1 
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Table 1. Characteristics of the clouds used to compare the models of dense molecular clouds described in the paper. The first column 
gives the name of the observed cloud, the second column the value of the visual extinction at the centre of the cloud, the third column 
the visual extinction around the cloud (extinction of the incident radiation field), the fourth column the lowest colour temperature 
measured in the cloud, the fifth column indicates whether the submillimetre emission of big grains is enhanced towards the cloud 
centre, the sixth column shows whether an anti-correlation between the emissivity spectral index and the colour temperature is 
observed, and the last column indicates the corresponding references. 



Object 



Central A w External A v 



Filament in Taurus < 10 ~ 0.4 - 0.5 

Filament in Musca ~ 8 ~ 0.3 - 0.4 

Planck -G126.6+24.5 - 3.5 <1 

Lynds 1780 4 ~ 0.25 

IC5146(core# 1) -21 \ 



12 1+"-^ 
- 12.5 
11.6±0.8 
14.9 + 0.4 
10.2 ±0.1 



(K) Submm [} -T References 



x3 4+"^ 


Yes 


\ 


Yes 


\ 


Yes 


- xl.5 


\ 


Yes 


\ 



Stepnik et al. 1 2003 ij jPlanck Collaboration| ( |2011d| 
Juvela et al. ( 201 1 ) 
Planck (Jollaboration1(l2 01 lb 
Kidderstad et al. ( ZOOoIr 



Kramer et al. I 



2003r 



Kidderstad & Juvela 



j20T0| 



Table 2. Dust model abundances and size distribution parame- 
ters (see Section 2.1 for details), p is the grain mass density, Y 
is the mass abundance per H, K2S{)^im is the opacity at 250 ;um, 
and p is the intrinsic opacity spectral index for A > \QQ fim for 
each dust population. For aggregates, the percentages give their 
degree of porosity (see Section[4!2]for details). 



ter is mostly distributed along elongated filaments in both the 
dense and the diffuse medium ( Abergel et al.|2010 Andre et al 



Small grains (DHGL) 




P 




ao 


Y 








(g/cm^) 




(nm) 




(cmVg) 




PAH 


2.24 


0.40 


0.64 


7.80x10-'' 


0.001 




SamC 


1.81 


0.35 


2.00 


1.65x10-'' 


0.002 


1.55 






Bi^ 


' grains (DHGL) 








P 


a 


ar,a, 


Y 


K250 iim 


P 








(nm) 








LamC 


1.81 


-2.8 


150.0 


1.45x10-3 


0.014 


1.55 


aSil 


3.5 


-3.4 


200.0 


7.8x10-' 


0.034 


2.11 






Accreted 


"rains 








P 


a 


ac,a, 


Y 




P 


accC 


1.81 


-2.8 


150.0 


1.6x10-3 


0.015 


1.55 


accSi 


3.22 


-3.4 


200.0 


8.6x10-3 


0.066 


1.52 


Aggregates 




P 


a 




Y 






0% 


2.87 


-2.4 


234.0 


1.02x10-- 


0.111 


1.33 


10% 


2.59 


-2.4 


242.0 


1.02x10-- 


0.140 


1.32 


25% 


2.16 


-2.4 


256.0 


1.02x10-2 


0.208 


1.30 


40% 


1.72 


-2.4 


276.0 


1.02x10-2 


0.331 


1.27 



201 01 lJuvela et al.|[20T0||201ll |Miville-Deschenes et al.||2010t 
Arzoum anian et al.|201 1). To account for this particular geome- 
try, we implemented a third type of cloud in CRT: finite circular 
cylinders. Our cylinders are divided into layers perpendicular to 
the axis of symmetry and each layer is divided into concentric 
rings. Each concentric ring is thus a cell of the modelled cloud, 
except for the central cell of each layer, which is a disk. One dis- 
crete cell has a constant density but density can vary from one 
cell to the other in both directions, that is to say as a function of 
the height along the axis of symmetry and as a function of the 
radial position. 

The CRT model consists of two independent parts. First, the 
radiation field is estimated at each position inside the cloud using 
Monte Carlo methods as described in |Juvela & Padoan| ( |2003[ l. In 
the case of cylinders, we consider incoming and outcoming pho- 
tons over the entire surface of the cloud, including the two caps. 
Second, using the radiation field and the selected dust model, 
dust temperature and emission are calculated. The CRT model 
allows us to add different populations of grains, with distinct 
spatial distributions. Dust temperature distributions are solved 
using DustEM (Compiegne et al.,201 1 ). The non-isotropic and 
multiple scattering events are taken into account, as well as the 
re-emission by dust grains by iterating CRT Puvela & Padoan 
.2003„Juvela.2005) . 



model abundances and the parameters of the size distributions 
are given in Tab. 2. We also consider the change in the grain op- 
tical properties towards dense molecular clouds, which are usu- 
ally assumed to be caused by their growth by accretion and/or 
aggregation. The populations of evolved dust are described in 
Section|4]and Tab. 2. 

2.2. Radiative transfer model 

We use the Monte Carlo rad iative transfer model described in 
Juvela & Padoanl ( |2003| and | Juvela] (2005 1, CRlf] (Continuum 
Radiative Transfer). Regarding the cloud geometry, until now, 
two types of clouds could be used to perform the calculations: 
spheres, divided into discrete cells with constant density, which 
are concentric shells, or three-dimensional clouds, divided into 
cubic cells that have a constant density. The second case is time- 
consuming as no assumption is made about the symmetries. 
Observations of the Galactic ISM indicate that interstellar mat- 



2.3. Definition of tine test cioud 

In CRT, a cloud is defined by its hydrogen density distribution, 
p(r). Cylindrical isothermal self-gravitating clouds ha ve distri- 
butions varying as p{r) cc in the outer regions (Ostriker 
1964). However, interstellar filamentary clouds usually exhibit 
much shallower profiles from p(r) oc r-'^ to r-^^, associated 
with a flat distribution at the centre ([Alves et al. 1998[ [Stepnik 



eraL][2003| [Arzoumanian et al.||20Tl| i. |Fiege & Pudritz] ( |2000l l 



demonstrated that these profiles are representative of magnetized 
filaments where the magnetic field influence dominates over 
gravity in the outer regions. For isothermal filaments, they found 
that p(r) oc r-'-^ to r^^ and for logatropic filaments p(r) oc r-' 
to r-'^. These shallower profiles may also be explained by non- 
isothermal, collapsing clouds, accreting the surrounding diffuse 
gas (Nakamura & Ume murajl 999) 1 . 

Dapp & B asu, ( ,2009) l and [Arzoumanian et al. ( 201 l| l showed 
that the column density profiles of interstellar filaments can be 
fitted with the following radial hydrogen density distribution 



3 Available at http ://wik i .helsinki . fi/display/-mjuvela@[ 
[heTsinki . f i/CRT| P^''> 



Pc 



1 +(r///o)-' 



(1) 
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where pc is the central density, which remains constant along the 
axis of symmetry, and Ho is the internal flat radius. We define 
the (pc, //o)-parameters for the clouds to be at equilibrium, with 
their ma s s per unit length equal to the critical value defined by 



Ostriker ( 1964j|, Mc,„ - 2c^/G, where G is the gravitational 
constant. Ihe sound speed, cs ~ 0.2 km/s, is assumed to be 
constant and is computed for a gas temperature of Tgas = 12 K. 
The steepness of the density profile, a, is usually in the range 
< 2.5 ( |Dapp & Basu 2009| [Arzoumanian et al.|20lT 



1.5 < a 

For our study, we assume that a - 2. The mass per unit length 
for the profile defined by Eq.[T]is 



M 



f 

Jo 



2np(r)rdr 



(2) 



M - npcHlXn 



1-Hl^ 

Ha 



(3) 



where R is the outer radius of the cloud. Assuming that M 
Merit, the (pc, ^^{))-parameters can be solved numerically 



npcHl In 



1 



2cl 



(4) 



We set ^? = 1 pc to be able to compare our results with those 
of Arzoumanian et al. (201 V\, who used R - 1.5 pc and those 
of 



Juvela et al.| (2012|, who chose R - 0.4 pc. The clouds 



are divided into 217 cells, consisting in 31 rings and 7 layers 



(see Section 2.2 1. We compute the (pc, //o)-parameters to model 
clouds with visual extinction at the centre from Ay - 1 to 20 
(Tab. 3). The emission maps produced have 41x41 pixels with 
a single pixel size of 0.025 pc, equivalent to 0.85' at a distance 
of 100 pc or 5" at 1 kpc. Finally, when not specified, these ide- 
alized filaments are illuminated by the isotropic interstellar stan- 
dard radiation field, ISRF, as defined by |Mafliis et"aL] ( [T983) l. We 
also study an alternative geometry for the clouds: Bonnor-Ebert 
spheres (Bonnor 1956 ; Ebert 1955^ . The results are presented 
in Appendix |B| and allow us to test the influence of the cloud 
geometry on the results. 



2.4. Dust colour temperature and emissivity spectral index 

Thermal emission from large dust grains dominates the observed 
emission in the far-IR and in the submillimetre. These grains are 
in thermal equilibrium with the interstellar radiation field and 
the observers often assume that their spectral energy distribution 
(SED) is well-approximated by a single grain size with a con- 
stant temperature and abundance along the line of sight 



)P colour 
ByiTcolour), 



(5) 



where 5v is the brightness. By is the Planck function, T colour 
is the dust colour temperature, eo is the dust emissivity at fre- 
quency vq (or the optical depth per unit column density tq/Nh 
in cm^/H), and PcoUmr is the emissivity spectral index. To com- 
pare our model with observations, we fit the model central 
pixel, which is also the brightest and the coldest pixel, with 
this modified blackbody where the factor Nh^o, the spectral in- 
dex jScoiour, and the colour temperature Tcoiour are free parame- 
ters. Three cases are considered. Except where otherwise stated, 
we either fit the entire SED from 100 jum to 3000 ^im (~ 100 
GHz), or perform the fit using the Planck-HFQ(High Frequency 

* We exclude the two HFI channels which are contaminated by CO 
emission, 100 and 217 GHz, as in Planck Collaboration 1 2011d i. 
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12 3 4 5 

A"' applied to ISRF 

Fig. 2. Colour temperature for the central brightest pixels (for 
SEDs fitted between 100 /im and 3000 jum) of the idealized fila- 
ments as a function of A"', which is the extinction applied to the 
ISRF (and the corresponding Go factor). The blue line shows the 
cylinder with Ai/ = 1 at the centre, green shows 5, red 10, light 
blue 15, and magenta 20. The true grain temperatures, Tdust, can 
be seen in Fig. [8] 



Instrument) filters at 857, 545, 353, and 143 GHz (350, 550, 
850, and 2 100 jum, respectively) , and the IRIsj^ 100 jum filter as 



( |2011d| , or we use the Herschel PACS0 



in Planck Coll aboration 

and SPIRrf] filters at 100, 160, 250, 350, and 500 ^lm. ExcepT 
in Section |3.2| where we consider the eff'ects of noise, we do a 
standard weighted least square fits, where all bands are given 
an equal weight, in order to determine TcoUmr and fi colour- The 
results do not depend on the absolute value adopted for the un- 
certainty (0.01% here), but only on the relative weighting of the 
data points. The values of Tcoiour and Pcoiour are also completely 
dependent on the wavelength range/instrumental filters used in 
the fitting procedure. This dependence can be seen in Fig.[T] for 
example, and comes from the different spectral distribution of 
the filters regarding the position of the peak wavelength of the 
big grains' emission. 

3. Dust with the properties of the diffuse high 
Galactic latitude ISM 

We first assume that in the ten dense idealized cylindrical fila- 
ments defined by Eq.[T]with the parameters listed in Tab. 3, dust 
has the same properties as in the diffuse ISM (Compiegne et al. 
201 1 ). The dust populations present in our modelled cylindrical 
clouds, with constant abundances, are thus PAHs, SamC, LamC, 
and aSil, described in Tab. 2 for DHGL. 



3.1. Effect of the radiative transfer 

As shown in Fig.[T];, a cloud with a visual extinction of Ay = 1 at 
the centre has a dust colour temperature higher than 17.5 K, irre- 
spective of whether the dust SED is fitted between 100 and 3 000 
fim or in the Planck HFI and IRIS 100 jum bands. When Ay = 10 
or 20, the colour temperature is lower and can be as low as either 



' IRIS is the new generation of IRAS (InfraRed Astronomical 
Satellite) images described in| Miville-Desch6nes & Lagache^ ( ^2005^ . 
* Photodetector Array Camera and Spectrometer. 
' Spectral and Photometric Imaging Receiver. 
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Table 3. Parameters defining the gas distribution of the modelled clouds. Central density (pc in H/cm^) and internal flat radius (Hq 
in pc) for the clouds with visual extinction at the centre Ay = 1 to 20, and for all the dust populations considered in the paper (see 
Tab. 2). 



DHGL (Dust at high galactic latitude) 


A„ 


1 


2 


3.5 


5 


7.5 10 


12.5 


15 


17.5 


20 


PC (H/cm3) 


460 


2900 


11480 


26 800 


69 300 135 350 


225 500 


341 000 


481500 


647 750 


Ho (pc) 


0.74 


0.16 


0.065 


0.039 


0.022 0.015 


0.011 


0.0090 


0.0075 


0.0063 


Accreted grains 


Pc 


370 


2 300 


9 350 


22000 


56 800 111 000 


185 600 


282 000 


398 000 


536 500 


Ho 


1.04 


0.18 


0.073 


0.044 


0.025 0.017 


0.013 


0.010 


0.0083 


0.0070 


Aggregates 0% 


Pc 


590 


3 620 


14100 


32 800 


84 500 164 500 


273 800 


413 000 


581800 


783 000 


Ho 


0.55 


0.13 


0.057 


0.034 


0.020 0.014 


0.010 


0.0081 


0.0067 


0.0057 


Aggregates 10% 


Pc 


520 


3 240 


12 700 


29 700 


76 500 149 000 


248 300 


375 000 


528 600 


711 100 


Ho 


0.64 


0.14 


0.061 


0.037 


0.021 0.014 


0.011 


0.0086 


0.0071 


0.0060 


Aggregates 25% 


Pc 


420 


2700 


10620 


24900 


64 350 125 800 


209 700 


317 500 


448 200 


603 400 


Ho 


0.84 


0.16 


0.068 


0.041 


0.023 0.016 


0.012 


0.0094 


0.0078 


0.0066 


Aggregates 40% 


PC 
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2200 


8 850 


20 800 


54 000 105 500 


176500 


268 000 


378 800 


510 800 


Ho 


1.24 


0.19 


0.076 


0.045 


0.026 0.018 


0.013 


0.010 


0.0085 


0.0072 



15.4 K or 14.6 K. In a cloud with high visual extinction, the ra- 
diation field at short wavelengths (visible/ultraviolet), which is 
mainly responsible for the grain heating, is more extinguished 
in the inner layers of the cloud. However, these colour temper- 
atures are still higher by several Kelvins than what is observed 
in dense molecular clouds (see Tab. 1). These dense clouds are 
often embedded in large HI or molecular complexes, leading 
to a lower and reddened incident radiation field at the surface 
of the cloud and thus to lower colour temperatures. Using the 
same clouds with central Ay = 1 to 20, we illuminate them 
with an extinguished isotropic ISRlj^ We extinguish the radi- 
ation by Af/' = to 5, corresponding to Go = 1 to 0.005, 
respectiveljp] The results, shown for fits to the dust SEDs be- 
tween 100 and 3 000 /urn, are presented in Fig. |2] A radiation 
field that is decreased by one magnitude of visual extinction, or 
with Go = 0.3, leads to a decrease in the observed colour tem- 
perature of less than 2 K. This diminution becomes smaller when 
the cloud central column density or visual extinction increases. 
Most of the UV/visible photons, which disappear when the radi- 
ation field is extinguished, have no influence on the inner layers 
of the densest clouds. Consequently, the extinction in the exter- 
nal radiation field alone cannot explain the colour temperatures 
measured in the objects described in Tab. 1 ( Bernard et al.|1992' 



|Stepnik et al.|2003] l. Similar results are obtained in the case of 
clouds model ed as Bonnor-Ebert spheres instead of cylinders 
(see Appendix [B). 

Fig.[T]also shows the Pcoiom — T colour relation obtained by fit- 
ting the dust SEDs between 100 and 3 000 /vm, and in the Planck 
HFI and IRIS 100 |/m bands. This relation exhibits a correla- 
tion instead of the observed anti-correlation between the emis- 
sivity spectral index and the colour temperature. Similar results 
were obtained by |Malinen et al.| ( 201 l| l, who modelled turbu- 
lent clouds containing gravitationally bound cores (see their Fig. 



As noted by 



Stamatellos & Whitworth 



1 2003 I, the radiation field 



illuminating embedded clouds is usually not isotropic but we keep this 
assumption for the sake of simplicity. 

' Go scales the radiation field intensity integrated between 6 and 13.6 
eV. The Mathis radiation field, ISRF, Go = 1, corresponds to an inten- 
sity of 1.6 X 10"^ erg/s/cm^. 



15a), and this correlation is also seen in the case of Bonnor-Ebert 
spheres (see Appendix [B). It can be explained by purely radia- 
tive transfer effects for grains with constant emissivity spectral 
indices. The clouds are optically thick at short wavelengths and 
the UV/visible photons, which heat the dust grains, are absorbed 
efficiently in the outer regions. Thus, the dust grains are colder 
at the centre of the clouds than on the surface. As the clouds 
are optically thin at long wavelengths, we observe the emission 
of all the grains present along the line-of-sight. The range of 
temperatures increases with the central density: the dust tem- 
perature at the edge of the clouds, which are illuminated by the 
ISRF, remains almost the same, regardless of the central den- 
sity, whereas it is significantly lower at the centre. Consequently, 
the resulting SED broadens as the extinction of the clouds in- 
creases, leading to a decrease in the emissivity spectral index to 
account for the flatter shape, while the colour temperature de- 
creases. This Pcoiour - Tcoiour Correlation, obtained using meth- 
ods including full radiative transfer calculations, contradicts the 
results of Shetty et al. ( 2009 1, who found that the line-of-sight 
temperature variations naturally lead to a Pcoiour — Tcoiour anti- 
correlation. This difference emphasizes the importance of the 
use of full radiative transfer calculations when modelling the 
dust emission from dense molecular clouds and is addressed in 
detail by |Juvela & Ysard| ( |2012| l. These authors demonstrated 
that the anti-correlation observed by 'Shetty et al. (2009 1 results 
from the particular temperature structure of the objects that they 
modelled, which are not likely to be representative of externally 
heated interstellar clouds. 



3.2. Effect of noise on the/3-T relation 

Radiative transfer effects lead to a correlation between the dust 
emissivity spectral index, PcoUmr, and the colour temperature, 
Tcoiour, whereas observations exhibit an anti-correlation (see Tab. 
1). [Shetty et"aL] ( [2009) and |Blain et"aL] ( [20031 ) showed that this 
anti-correlation can be produced by intrinsic noise in the obser- 
vations, which produces a degeneracy between Tcoiour and PcoUmr- 
This anti-correlation has been observed in various interstellar 
environments by Planck and Herschel instruments, for instance 
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A T , (K) 

V colour 

Fig. 1. Case of DHGL populations, a) 100 fim brightness map of the cloud with Ay - 10 at the centre, b) The black solid line is the 
SED of the central pixel and the black circles are the same SED integrated in Planck-HFI bands (143, 353, 545, 857 GHz) and IRIS 
100 yum band. The green dashed line is the best fit between 100 fim and 3000 fj.m, and the green diamonds show the best fit for the 
SED integrated in the IRIS 100 /im and in the Planck HFI bands, c) Colour temperature as a function of visual extinction for the 
central pixel for a series of models with different central Ay (see Tab. 3). The blue line shows the colour temperatures obtained by 
fitting the SEDs between 100 fim and 3000 fim, and the green line when the fit is made in the Planck-HFI and IRIS bands as defined 
in section[3] d) Emissivity spectral index, /3a,iour, as a function of the colour temperature. Line styles are the same as in the previous 
figure. y6(LamC) and y6(aSil) are the intrinsic opacity spectral indices of the LamC and aSil populations. The true grain temperatures, 
Tdusi, can be seen in Fig.|8] 



by|Planck Collabor ation (201 ld'c\'Paradis et al. (7 0T0| |Rod6n 
let al.| ( |2010| l, |Juvela et al. (201 1 ) and Bracco et al. (20Tl] ). These 
authors suggest that the observed Pcoiour - T colour anti-correlation 
could be produced mostly by intrinsic variations in the dust op- 
tical properties with temperature rather than noise. 

To estimate qualitatively the extent to which noise could lead 
to an apparent PcoUmr - Tcoiom- anti-correlation, we examine the 
effect of observational noise on the {Pcoiom—T colour) relation with 
the help of Monte Carlo simulations. We consider the ten clouds 
described in Tab. 3, assuming that they are all located at a dis- 
tance of 100 pc, with DHGL grain populations, and we use the 
associated surface brightness maps for diff'erent external radia- 
tion fields: ISRF and six values of attenuation of the external 
radiation field corresponding to the ISRF extinguished by A"' 
up to 5 (Fig.[2]i. 

We compute 1 000 realisations of the observations of the 
modelled clouds. In each realisation, the surface brightness val- 
ues are modified by the addition of normally distributed obser- 
vational noise. We consider the case of observations obtained 
with the combination of Planck-HFI and IRAS 100 fim bands. 



for which the modelled brightness maps are smoothed to a com- 
mon resolution of 5'. We also study the case of observations 
made with the Herschel PACS and SPIRE instruments for a com- 
mon resolution of 37". The standard deviation in the noise is 
assumed to be 0.06, 0.12, 0.12, and 0.08 MJy/sr in the IRAS 
100 fim band and in the Planck HFI bands at 857, 545, and 353 
GHz, respectively ( Planck Collaboration 201 la ). In the case of 
Herschel, absolute values of noise equal to 8.1, 3.7, 1.2, 0.85, 
and 0.35 MJy/sr (at the beam resolution) are used for the five 
bands at 10 0, 160, 250, 350, and 500 yum, respectively ( [Malinen 
|et al. [20lT| . For each realisation (observations of all frequen- 
cies, including the noise), the Tcoiour and PcoUmr values are esti- 
mated using a standard method for the central pixels of the 
series of filaments described in Tab. 3. Their uncertainties come 
from Monte Carlo simulations. These results are used to simulate 
surveys that would equally cover clouds from all the examined 
cloud categories (with Ay = 1 to 20 and variable attenuation of 
the ISRF A"' = to 5). For each modelled cloud, we extract 
200 random realisations, i.e., 200 Tcoiour and Pcoiour values de- 
rived from observations of that model with 200 different Monte 
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10 15 20 

^colour 



25 



300 




0.55 



0.60 



0.65 



Fig. 3. a) and b) Two examples of the (Tcohur, /^colour) distribu- 
tions for the observations simulated in the Planck-HFI and IRAS 
100 yum bands for the central pixels of the series of filaments de- 
scribed in Tab. 3. The red squares show the results of the fits 
for the clouds described in Section 3.1 containing DHGL grain 



populations and illuminated by the ISRF, extinguished or not, 
without noise (10 clouds x 7 radiation fields). The black crosses 
show the same after adding the estimated statistical noise to the 
simulated data (the length of the arms of the black crosses are the 
error bars). The green dashed lines are the power-law fits of these 
two particular realisations and the green solid line corresponds 
to the average relation over the 1 000 cases. For clarity, t he p lot 

for 



3.2 



includes only every 50* of the 14 000 points (see Section 
details). The black line is the fit to the observed Pcoiour - Tcoiour 
relation measured by |Planck Collaboration ( 2011c I towards cold 
clumps with Planck HFI and IRAS 100 yum bands, c) Probability 
distribution of the B parameter obtained from 1 000 cloud sam- 
ples, similar to the two examples in th e up per frames (with B 
defined as ficoiour - , . see Section 3.2 for details). 



Carlo estimates of the observational noise. This leads to a sam- 
ple of 14000 {Tcoiour, Pcoiour) points (10 clouds X 7 radiation 
field intensities x 200 noise realisations). Figures [3^ and b, and 
Figs. |4^ and b show four examples of these realisations, where 
for clarity only every 50''' of the 14000 points are plotted. The 
red dots show the results obtained when noise is not taken into 
account: for comparison, the green line of Fig. [TJl, obtained for 
the ten modelled clouds illuminated by the ISRF, corresponds 
to the red dots that are at the extreme right of Figs. [3^ and b. 




350 



10 15 20 

^colour f"*^' 



25 




0.9 -0.6 -0.3 0.0 0.3 0.6 0.9 1.2 
B 

Fig. 4. a), b), and c) Same as Fig. [i] for the observations simu- 
lated in the Herschel bands for the central pixels of the series of 
filaments described in Tab. 3. In a) and b), the black line is the 



fit to the observed Pcoiour — Tcoiour relation observed by |Paradis 
et al. (2010 1 in the Galactic plane with Herschel. 



Then moving to the left, the red dots correspond to decreasing 
intensities of the radiation field or similarly to increasing A"'. 
We then proceed to fit the 14000 {Tcoiour, ^colour) points with a 
relation PcoUmr - ^T^^^/o,,, ' which is commonly used to fit real 
observations, using a standard method, where the points are 
weighted according to their error estimate^^ The sign of the B 
parameter indicates whether the Pcoiour - Tcoiour relation is cor- 
related {B > 0) or anti-correlated {B < 0). These procedures, 
the random selection of the 14000 points and the fitting of the 
power-law, are repeated 1 000 times to get information about the 
distributions of the A and B parameters (Figs.|3j; and|4j;). 

In the case of Planck (h-IRAS), Fig.[^ shows that 0.62 < B < 
0.68. Even if flattened, the relation between Pcoiour and Tcoiour 
remains positively correlated with (^colour ^ T^ohur- '■^^ ^^^^ 
of Herschel, the values of B are more scattered. However, steep 
relations with B < -1 are excluded at a 99.9% confidence level, 
while relations steeper than B - -0.5 are found for only ~ 10% 
of the cloud samples. The difference between the Planck and 



In Figs.jsji and b, and Figs. [4^ and b, the outliers with low colour 
temperatures and high emissivity spectral indices are not captured by 
the fits (green lines). These points indeed represent fewer than 3% of 
the total number of points in the sample in the case of Planck, and fewer 
than 5% in the case of Herschel, and have large error bars. 
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Herschel cases comes from the high level of noise introduced in 

the PACS channels^ 

[Planck Collaboration|P01 1^ used a combination of Planck- 
HFI data (857, 545, and 353 GHz) and IRIS 100 ^m^ data to 
determine the dust properties in the starless cold clumps de- 
scribed in the Cold Core Catalogue of Planck Objects (C3PO). 
They observed a n anti-cor relation that could be fitted with 
Pcoiour ~ T colour- Pl™ck Collaboration (2011d) found a simi- 
lar Pcoiour - Tcoiour relation m the launis molecular complex, 
using Planck-HFI and IRAS bands. In the case of Herschel, 
[Paradis et al.| ( |2010| l measured an anti-correlation that c ould be 
fitted with Pcoiour ~ 7^co';,^„. in the Galactic plane, while Bracco| 
et al. ( 201 1\ found a comparable relation in a Galactic cirrus. 



10" 



Consequently, if we assume that our sample of modelled clouds 
is representative of the ranges of dust temperature and opacity 
spectral index of the observations, noise as the unique expla- 
nation of the observed anti-correlation appears to be unlikely. To 
test the influence of the cloud geometry, we did the same analysis 
using externally heated Bonnor-Ebert spheres and found similar 

results (see Appen dix [B]). 

As pointed by |Shetty et"aL] ( [2009| and |Blain et"aL] ( |2003) , 
intrinsic noise in the data produces a degeneracy between Pcoiom 
and Tcoiour potentially leading to an anti-correlation (Fig.|4]l or at 
least a flattening of the relation (Fig. [3]). Nevertheless, in our sim- 
ulations noise alone cannot explain the strong anti-correlation 
observed. For starless interstellar clouds, it seems probable that 
the anti-correlation is explained instead by variations in the dust 
optical properties with some contribution from the noise (see 
Section [S!2j l. However, if the observers do not exclude sources 
with very low signal-to-noise ratios (S/N's) (e.g., S IN < 4 at 
250 jum, see Appendix B), the noise can then have a dominant 
role on the observed Pcoiom - Tcoiour relation. 



4. Evolution of grains with medium density: 
accretion and aggregation 

The size of the dust grains is expected to increase with the 
density of the medi um (|Ossenkopf||1993 [Bernard et al.||1999 



[Stepnik et al.[|2003| [Steinacker et al.|[20i6| l, hence we investi 
gate the efifects of grain growth on the long-wavelength SED 
of dense clouds. The parameters we consider when comparing 
with observations are the emission deficit in the mid-IR, the 
higher emissivity in the submillimetre, the decrease in the grain 
colour temperature, and the Pcoiour - Tcoiour relation. The deficit 
in the mid-IR is usually explained by a lower abundance of small 
stochastically heated grains (PAH and SamC populations), and 
the higher submillimetre emissivity by coagulation ( Ossenkopf 



& Henning 1994 Stognienko et al. 1995 i. Four scenarii for grain 



growth can be considered: 

- aggregation of small grains with small grains ; 

- aggregation of big grains with big grains ; 

- accretion of small grains onto big grains (hereafter accretion 
scenario) ; 

- aggregation of all grains with all grains (hereafter aggrega- 
tion scenario). 

The first scenario, aggregation of small grains with small grains, 
allows us to reproduce the emission deficit in the mid-IR but 
has no significant impact on the grain submillimetre emissivity: 
large grain submillimetre emission and colour temperature are 
almost unchanged in this scenario. The second scenario leads to 
a decrease in the big grain colour temperature and an increase in 
the submillimetre emissivity, but does not explain the emission 
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Fig. 5. Absorption efficiencies, Qats, for an aSil grain (dashed 
line), a LamC grain (dotted line), and an accSi grain (solid line), 
with sizes as defined in section [4T| 



deficit in the mid-IR. The two last scenarii both explain qual- 
itatively the emission deficit in the mid-IR with the loss of all 
the small grains. We now describe in detail their eff'ects on the 
dense cloud SEDs. To perform this analysis, we consider that 
the ten cylindrical clouds described in Tab. 3 contain only either 
accreted grains (Section 4. 1 1 or aggregates (Section 4.2 1. 



4.1. Accretion 

All the small carbonaceaous particles (PAHs and SamC) are ac- 
creted at the surface of larger grains (aSil and LamC) as uniform 
carbon mantles of the same thickness, regardless of the size con- 
sidered, to form accreted carbonaceous particles and accreted 
silicates (accC and accSi, see Tab. 2.) This increases the total 
mass of large grains by 10% and their volume by 17%, while 
their radius remains almost unchanged ( — 1-6%). 

The LamC and accC populations have the same optical prop- 
erties and size distribution. In contrast, accreted silicates, accSi, 
are composed of 84% of aSil and 16% of carbon in volume (the 
amC mantle has a specific density of 1.81 gjcvD?, compared with 
3.5 glcvc? for silicates, see Tab. 2). Their optical properties are 
calculated with the effective medium theory (EMT), using the 
Bruggeman rule ( Bohren & Huffman 1983 ). The dielectric func- 
tion of the composite medium, egff , is the solution to Eq. 6 involv- 
ing the dielectric function of the first component, e\, of fractional 
volume /, and the second component, €2, of fractional volume 
1 - / given by 



/ 



feff 



+ (!-/) 



■ ^eff 



62 + 2eeff 



0. 



(6) 



For a two-component mix similar to that considered here, sili- 
cateous and carbonaceous, this equation is simply quadratic in 

The variation in the optical properties of the accreted grains 
compared with aSil and LamC are shown in Fig. |5] for grains 
containing the same amount of silicate (oaSii = 0.1 ;um o 
flaccSi = 0.106 yum). The absorption efficiency of accSi grains, 
Qabs, in the far-IR and the submillimetre is flatter {[i - 1.52) 
than for aSil (j6 = 2.11) and LamC (/? = 1.55): at 250 jum, it is 
multiplied by a factor of 1 .54. Qats also increases in the visible 
and in the near-IR. As a result, accSi grains absorb more energy 
than aSil grains. Both their submillimetre emissivity and their 
temperature are higher than for DHGL dust (Fig. |7^). Despite 
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of composite aggregates, heavy computational methods such as 
PDA or T-MA TRX must be used (Stognie nko et al.|1995[|Fogel 
& Leung|1998; Mackowski 2006 , Kohler et al.|201 1 ' 



Wavelength (|jni) 

Fig. 6. Absorption efficiencies, Qats, for an aSil grain (black 
dashed line), for a LamC grain (dotted black line), for an ag- 
gregate without porosity (blue), with a porosity degree of 10% 
(red), 25% (green), and 40% (magenta), with radii as defined in 
section |4!2] All these grains have the same volume of material, 
excluding the void only responsible for the increase of the size. 



being enhanced, their emissivity at long wavelengths is indeed 
not high enough to counterbalance the larger amount of energy 
absorbed at shorter wavelengths. This increase in the colour tem- 
perature does not match the observations. We also note that the 
accSi grains have a lower spectral index value that is, however, 
correlated with the colour temperature (Fig.|7]3). The results are 
qualitatively similar for Bonnor-Ebert spheres (see Appendix[B]i. 

4.2. Aggregation with or witliout porosity 

In this scenario, all grains are aggregated together to form 
mixed aggregates of small and big grains, which are necessarily 
porous. These aggregates are expected to form in dense molec- 
ular clouds, possibly with large fractions of void or porosity 
( [Dorschner & Henning|[T995| ). In our simple model, the frac- 
tion of silicate, carbon, and vacuum in the grains do not de- 
pend on their sizes. It is defined using the fractional porosity 
parameter p. The aggregates are composed of 63(1 - p)% sili- 
cate, 37(1 - p)% carbon, and 100/:'% vacuum. These proportions 
correspond to the conservation of the mass during the aggrega- 
tion of the DHGL grain populations defined in Tab. 2. As the 
total volume of big grains is dominated by silicates, the size 
distribution is modified starting from the aSil size distribution. 
Grains are on average larger and the size distribution is multi- 
plied by a factor (0.63 x (1 - p)y^^^. From a^^-'^, the power-law 
then becomes a"^ '* to mimic the loss of the small grains by ag- 
gregation, which allows us to reproduce the trend of an increase 
in Ry towards dense media (Fitzpatrick & Massa 1988; Mathis 
& Whiffen|| 19891 ICardelh & craytonp"991| [Campeggio et'aLl 

The optical properties of this new grain population are com- 
puted with the Bruggeman rule, as presented in Section 4. 1 but 
with three components of silicate, carbon, and void (e = 1). We 
first mix the carbon and silicate materials as in Section [4~Tl with 
their respective volume fractions of 63% and 37%. We then mix 
the resulting grains with void assuming the same two-component 
Bruggeman rule. This theory of an effective medium does repro- 
duce the general trend that we want to test here, where the inclu- 
sion of porosity is a way to increase the submillimetre emissiv- 
ity of grains. For a precise calculation of the optical properties 



The variation in the optical properties of the different popula- 
tions of aggregates compared with those for the aSil and LamC 
grains are shown in Fig. |6] We compare grains containing the 
same volume of material, meaning that the void is excluded from 
the calculation of the volume for the aggregates. The aggrega- 
tion of carbonaceous and silicate grains increases the absorption 
efficiency in the near-IR and the visible with respect to pure sil- 
icate cores. The trend is similar in the submillimetre: the sub- 
millimetre emissivity is increased first by the mixing between 
carbons and silicates, and second by the addition of porosity, 
a well-known effect ( |Voshchinnikov et al.||2006| . The increase 
in the emissivity for the aggregates wins over the increase in 
absorptivity, so that the colour temperature of the aggregates is 
lower and decreases with increasing porosity (Fig. |7^). In addi- 
tion, the emissivity spectral index in the submillimetre decreases 
from j6 = 2.11 for aSil, to 1.33 for aggregates without poros- 
ity, 1.32 with 10% porosity, 1.30 with 25% porosity, and 1.27 
with 40% porosity. Through aggregation, the dust emissivity is 
increased by a factor equal to the ratio of the Qabs of the ag- 
gregates to the mean Qabs of isolated DHGL carbon and silicate 
grains of the same volume. As carbon and silicate 0.1 fim grains 
have almost the same Q„bs at 250 fim (Fig. |5]l, the increase in 
the dust emissivity brought about by aggregation is simply the 
ratio of the Qafe.s of the aggregates of Fig. |6]to the Q„i,s of a 0.1 
pm silicate (or carbon) grain. We find an increase of 1.84 for 
a grain without porosity, 2.16 for a grain with 10% of poros- 
ity, 2.84 for 25%, and 3.89 for 40%. When the modification of 
the size distribution is taken into account, these factors are even 
higher: in the case of DHGL populations, only 83% of the to- 
tal volume of grains emit in the submillimetre (excluding PAHs 
and SamC), whereas 100% do in the case of aggregates (PAHs 
and SamC are now aggregated with big grains). Thus, the factors 
become 2.2, 2.5, 3.3, and 4.6 for aggregates with a porosity frac- 
tion of 0, 10, 25, and 40% respectively. For a cloud with a central 
Ay - 10, this leads to a higher "observed" emissivity by a factor 
of 3.1, 3.6, 4.8, and 6.9, for grains with a porosity degree of 0, 
10, 25, and 40% respectively, when the dust SEDs are fitted in 
Planck bands. These values agree well with what is observed in 
dense molecular clouds (see Tab. 1). We also note that the "ob- 
served" emissivity spectral index, ficoiour, is again correlated with 
the colour temperature in an opposite way to what is observed in 
the ISM (Fig.|7]5). Similar results are obtained for Bonnor-Ebert 
spheres (see Appendix [B). 

The aggregation of grains in mixed aggregates of small 
and big grains is able to explain simultaneously the deficit of 
emission in the mid-IR, the low colour temperature, and the 
higher submillimetre emissivity observed in dense molecular 
clouds (see Tab. 1). However, grain growth cannot reproduce the 
Pcoiour - T colour anti-corrclation, which might originate in grain 
intrinsic properties and in particular in variations in the emis- 
sivity spectral index with the wavelength and/or the temperature 
P{A, T) (see Section|5]|. 

4.3. Colour temperature versus grain temperature 

In dense molecular clouds, the intensity of the radiation field de- 
creases from the edge to the centre. Consequently, dust grains 
have different temperatures for different positions inside the 
cloud and these temperatures differ from the measured colour 
temperature. Figure[8]shows the relation between the colour tem- 
perature, Tcoiour, and the temperature of the grains, Tdusi, which 
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colour ' Radius (pc) 

Fig. 8. a) Equilibrium temperature of dust in the central cell of the modelled clouds, Tdust, averaged over the grain size distribution, 
as a function of the colour temperature, T colours calculated by fitting the SEDs between 100 yum and 3000 ^m. The DHGL case is 
represented by the black lines (the solid line is aSil and the dashed line is LamC), the accretion case by the cyan lines (the solid line 
is accSi and the dashed line is accC), the coagulation case with 0% of porosity by the blue line, with 10% of porosity by the red line, 
with 25% by the green line, and with 40% by the magenta line, b) As in Fig. [8^) but at the edge of the cloud, c) Tdust as a function 
of the cloud radius for a cloud with central Ay - 10. The colour code is the same as in|8^). d) The same for the DHGL case (the 
circles show aSil and the diamonds LamC) for a cloud illuminated by the standard ISRF (dashed line) and illuminated by the ISRF 
extinguished by one visual magnitude (solid line). 



is defined as the individual equilibrium temperatures averaged 
over the grain size distribution. We determine one dust tempera- 
ture per cell and per grain population. 

Fig. presents the dust temperature at the centre of the 
clouds as a function of the colour temperature obtained by fit- 
ting the dust SEDs between 100 and 3 000 /vm. Except for the 
purely carbonaceous grains (LamC), T^ust is always lower than 
Tcoiour for clouds with Ay > 1- This is the result of the mixing 
of all the dust temperatures from the warmer edge to the colder 
centre of the clouds. The relation between the central r^,,,, and 
Tcoiour is almost linear for the clouds with central Ay < 10. 
This means that for this kind of moderately dense clouds we 
should be able to estimate the physical dust temperature at the 
centre, assuming one dust population and knowing the colour 
temperature. The results of these linear fits are given in Tab. 4. 
However, we expect these relations to be more complicated if 
several dust populations are mixed along the line-of-sight, for 
instance DHGL grains at the edge and aggregates at the cen- 
tre if we assume that the grains evolve within dense molecular 
clouds. FigurejS}?) shows the dust temperature at the edge of the 
clouds as a function of the colour temperature. For Ay = 1 to 20, 
Tdusi does not vary by more than 0.5 K and, as expected, its val- 
ues remain close to the dust temperature measured for the same 
population of grains illuminated by the ISRF, without radiative 
transfer However, we note a slight increase in Tdusi when Tcoiour 
decreases, or similarly when the central Ay increases. For varia- 



tions in Tcoiour from about 2.5 K to 3.5 K, depending on the dust 
population considered, Tdusr increases by ~ 0.5 K. This surpris- 
ing rise in the dust temperature can be explained by the decrease 
in the central flat radius, //q, with increasing central density, pc 
leading to a slightly higher incident radiation field in the outer 
layers (Tab. 3). 

Fig. [8]: displays the dust temperature as a function of the 
radial position in the cloud, for a cloud with a visual extinc- 
tion at the centre equal to 10. The temperature depends on the 
grain composition and it decreases from 7 K to 9.5 K from the 
edge to centre, in agreement with previous studies ( Bernard et al. 
[19921 [Evans et al.||200Tl [Zuccoiii et ai:i|200Tl |Fischera||201 iT 
For the central layers of the cloud, the dust temperature is al- 
most constant and equal to from ~ 8.5 K to 14 K depending on 
the dust population, the hottest grains being the carbonaceous 
ones (LamC). This is explained by grains being mostly heated 
by the visible and UV photons of the incident radiation field, 
which are mainly absorbed in the outer layers of the cloud. This 
behaviour is found for all the clouds we modelled and is the 
most visible when the central Ay is the highest. We also note 
that the decrease in the dust temperature from the edge to the 
centre of the cloud is smaller when the incident radiation field is 



extinguished (Fig. |8p), as similarly observed by Stamatellos & 
|Whitworth| ( |2003| ) and |Fischera| ( |201 l| l. The more the radiation 
field is extinguished, the smaller the temperature variation. 
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Table 4. Linear fits between the dust temperature at the centre 
of the clouds as a function of their colour temperature and their 
dust content, for clouds with a central Ay < 10. 
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Fig. 7. a) Colour temperature as a function of visual extinction 
for the central pixel. Results are shown for accreted and aggre- 
gated grains as defined in Tab. 2. Solid lines show the colour 
temperature calculated by fitting the SEDs between 100 yum and 
3000 fim, and the dashed lines in the Planck-HFI bands as de- 
fined in section 12.41 The black line shows the case of accreted 
grains, the blue line of aggregates without porosity, red of a 
porosity degree of 10%, green 25%, and magenta 40%, as de- 
fined in Tab. 2. The true grain temperatures, Td„si, can be seen 
in Fig. [8] b) Emissivity spectral index as a function of colour 
temperature. Line styles are the same as in the top figure. The 
numbers in parenthesis are the intrinsic opacity spectral indices 
of the different grain populations. 



5. Dust intrinsic properties 

Grain growth and purely radiative transfer effects can explain 
neither the Pcoiour — T colour anti-correlation, nor the variations in 
the emissivity spectral index with wavelength that are observed 
in the ISM. These two observations may thus originate in varia- 
tions in the grain intrinsic properties. Hence, we study the effects 
of variations in the optical properties of the grains with wave- 
length, /?(/!), and then temperature, y6(r^„j,), on the dust SEDs of 
the clouds. 



5.1. Variations in optical properties with waveiengtti: /3(A) 

We investigate the impact of spectral variations in the grain op- 
tical properties on the emerging emission from dense molecular 
clouds. We perform this analysis using laboratory data in the 
FIR and the submillimetre to modify the absorption efficiencies, 
Qabs, of LamC and aSil grains. We consider samples of amor- 
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Fig. 9. a) Absorption efficiencies for LamC (solid line), BE 
(dashed line), and ACAR (dotted line), for grains with a = 0.1 
fim. The small window shows the ratios of the absorption effi- 
ciencies of BE and ACAR to LamC with the same linestyles. 
b) Absorption efficiencies for aSil (solid line), ENSTA (dashed 
line), and OLIVA (dotted line), for grains with = 0.1 //m. The 
small window shows the ratios of the absorption efficiencies of 
ENSTA and OLIVA to aSil with the same linestyles. 



phous materials, which are representative of interstellar grains 
( |Draine|2003a[ [Kemper et al.|2004[|Fitzpatrick & Massa|2007 1. 



Mennella et al. ( 1998) measured the absorption coefficient of 
two different amorphous carbonaceous cosmic analogue grains 
between 20 yum and 2 mm at 24 K. Their first sample was 
composed of a population of chainlike aggregates of spheroidal 
grains with a diameter around 10 run (ACAR sample), and did 
not contain hydrogen. The second sample consisted of a soot 
with the same morphology but for spheroidal grains with a di- 
ameter of around 30 nm (BE sample). The BE sample is repre- 
sentative of ip^-rich hydrogenated amorphous carbon. To build 
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10" 
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Fig. 10. a) SEDs for the central pixel of the cloud with central 
Ay = 10 for a mixture of LamC and aSil (black), BE and ENSTA 
(green), BE and OLIVA (magenta), AGAR and ENSTA (blue), 
and AGAR and OLIVA (red), b) Emissivity in arbitrary units for 
the cloud with central Ay = 10, for a mixture of AGAR and 
OLIVA (solid grey line), and a mixture of AGAR and ENSTA 
(solid black line). Absorption efficiencies in arbitrary units for 
grains with a diameter of 0. 1 fim: AGAR (dotted line), ENSTA 
(dashed line), and OLIVA (dashed-dotted line). 



the absorption cross-sections of the two dust populations, for 
wavelengths shorter than 20 /im, we use the absorption effi- 
ciency of the LamG population and normalise the AGAR and 
BE absorption coefficients at 20 fim for longer wavelengths. For 
wavelengths longer than 2 mm, we then extrapolate Qahs using a 
unique spectral index as measured at 2 mm: /3(2 mm) - 0.95 for 
the AGAR population and 1.30 for the BE population. The re- 
sults are shown in Fig.|9^). Both the AGAR and BE populations 
have a flatter Qahs than the LamG population in the far-lR and in 
the submillimetre (for A > 20 //m). At 250 fim, when compared 
with LamG, the absorption efficiency is higher by a factor of 2.2 
for an AGAR grain with a = 0.1 fim and by 2.9 for a BE grain 
with the same size. The BE grains have a higher absorption ef- 
ficiency ffian the AGAR grains for wavelengths between 20 and 
3000 fim. For the AGAR and BE populations, we use the same 
size distribution and abundance as used for LamG grains. 

For the silicates, we use the data from [Goupeaud et al.| 
( 201 l| l, who measured the opacity of interstellar dust analogues 
at 10 K between 100 fim and 1500 fim. Two of their sam- 
ples are considered. The first is composed of amorphous sili- 
cate with a stoichiometry close to olivine (Mg2.3Si04, sample 
Fl in Goupeaud's nomenclature, OLIVA hereafter) and the sec- 
ond of amorphous silicate with a stoichiometry close to enstatite 



(MgSi03, sample E, hereafter ENSTA). In both cases, for wave- 
lengths shorter than 100 fim, we use the absorption efficiency of 
aSil and for longer wavelengths, we normalise the OLIVA and 
ENSTA absorption coefficients at 100 yum to produce the absorp- 
tion efficiencies. Finally, for wavelengths longer than 1 500 /im, 
we extrapolate Qabs using a unique spectral index as measured at 
1 500 fim: /3(1 500 fim) = 4.30 for OLIVA and 0.95 for ENSTA. 
The results are shown in Fig.|9]3. The emissivity spectral index of 
the absorption efficiency of OLIVA changes from yS = 2.1 below 
~ 700 fim to y6 = 3.6 above, and the spectral index of ENSTA 
from /3 = 2.5 for A < 600 fim, to yS = 1.7 for 600 < ^ < 800 fim, 
and to p - 0.9 for longer wavelengths (Goupeaud et al. 2011, 
see their Tab. 3). These indices deviate from the aSil emissiv- 
ity spectral index, which is equal to 2.11 for A > 100 fim (Fig. 

For wavelengths longer than 1 200/1 500 fim, the absorp- 
tion efficiency of OLIVA/ENSTA is extrapolated and the dif- 
ference from aSil, although very pronounced, is questionable. 
However, this wavelength range does not influence the results of 
the models (see the following paragraph). 

We now consider that the ten cylindrical clouds described 
in Tab. 3 contain various mixtures of the AGAR, BE, OLIVA, 
and ENSTA grains (for more details, see also Appendix [A). The 
resulting dust SEDs for these different mixtures of grain popula- 



tions are shown in Fig.[TO}i (AGAR or BE + OLIVA or ENSTA). 
The colour temperatures and emissivity spectral indices are de- 
termined by fitting the dust SEDs between 100 and 3 000 fim 
(Fig.[TT]i. The four mixtures produce lower colour temperatures 
than the DHGL populations, which is due to their increased 
emissivity in the far IR. For a central Ay - 10, the colour tem- 
perature is decreased by ~ 0.2 K for a mixture of AGAR and 
ENSTA, ~ 1.2 K for AGAR and OLIVA, ~ 1 K for BE and 
ENSTA, and ~ 1.9 K for BE and OLIVA. The temperature de- 
crease is the strongest for the dust populations which have the 
highest absorption efficiencies in the far-IR and the lowest at 
longer wavelengths {A > I 000 - 2 000 fim). For the dust temper- 
ature, the important spectral range is thus the far-IR, where large 
grains emit most of their energy in dense molecular clouds. This 
means that the extrapolation of the grain optical properties above 
2 mm for the carbonaceous grains and 1 500 fim for the silicates 
has no influence on the estimate of the dust colour temperature. 
However, Fig. 1 1 shows that the different values of /3 in the sub- 



millimetre for OLIVA and ENSTA lead to a flatter spectral index 
for mixtures containing ENSTA (1.0 < yS < 1.2) than for OLIVA 
(1.2 </?< 1.4). 

For data analysis, it is useful to know whether the nature 
of the dust grains can be determined by observing the spectral 
variations in the SEDs of the clouds. To address this question, 
we calculate the dust emissivity using the colour temperatures 
estimated previously (Fig. 1 1 1 and the hydrogen column density 
Nh 



BAT,oiour)NH Nh 



in cm /H. 



(7) 



Fig. 10 shows the example of mixtures of AGAR with either 
ENSTA or OLIVA grains for a cloud with Ay - I. The results 
are highly dependent on the grain populations included in the 
model. In the case where the silicates are assumed to be OLIVA- 
like, the slope of the emissivity of the resulting SED is similar 
to that of the carbonaceous grains for A > 300 - 1 200 fim de- 
pending on the cloud density, and regardless of the nature of 
the carbonaceous grain population (AGAR, BE, or LamG). This 
can be explained by the strong decrease in the emissivity of the 
OLIVA grains relative to that of the carbonaceous grains at long 
wavelengths (Fig. |9]l. Shorter wavelengths are then too close to 
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Fig. 11. Emissivity spectral index as a function of colour tem- 
perature for fits of the dust SEDs between 100 and 3 000 jim for 
the central pixels of the series of filaments described in Tab. 3. 
Line styles are the same as in Fig. 10 1). The numbers in paren- 



thesis are the intrinsic opacity spectral indices of the different 
grain populations. 



the peak of the SED to be used (Fig. 10 1)- In this spectral range 
(A < 300 i-im), the spectral shape of the SED is strongly influ- 
enced by radiative transfer effects (i.e. by temperature mixing) 
and consequently no spectral variations can be associated with 
the grain optical properties. When the silicates are assumed to 
be ENSTA-like, the results are then less clear. First, for a mix- 
ture with BE grains, we cannot make any obvious correlation be- 
tween the slope of the absorption efficiencies and the SED emis- 
sivity, except for 800 < A < I 700 /im. In this spectral range, 
the slope is similar to the one of ENSTA but differs at longer 
wavelength owing to the strong handover of the BE absorption 
efliciency around 2 000 ;um. Second, for a mixture with ACAR 
grains that have the same emissivity spectral index as ENSTA 
grains for /I > 2 000 yum, the emissivity spectral index is simi- 
lar to the one of ACAR for A > 720 fim. It is also possible to 
distinguish the bump of the ACAR absorption efficiency around 
1 400 fim. We conclude that it could be possible but quite diffi- 
cult to deduce the nature of the dust grain populations from the 
spectral variations in the dust emissivity estimated from the dust 
SED. As seen in Section 3.2 the observed dust emissivity spec- 
tral index may also be affected by noise, meaning that even more 
confusion should be expected in real observations. 

We also note that variations in the emissivity spectral index 
with wavelength alone can explain neither the strong decrease in 
the colour temperature towards dense molecular clouds nor the 
/^colour - T colour anti -Correlation (see Fig. 111. The last one may 
be caused by intrinsic variations in the grain optical properties 
with temperature. 



5.2. Variations in optical properties with temperature: piT) 

To investigate the influence of intrinsic variations in the emis- 
sivity spectral index with the temperature, /ig^^X^dust), we use 
the temperature-dependent grain emissivity tool of DustEM de- 
scribed by Compiegne et al. (2011), which allows us to vary 
the spectral index of the absorption efficiency, ySg^,^,, with the 
grain temperature, Td„st- In DustEM, the temperature depen- 



Fig. 12. The solid line is the Pcoiour — Tcoiour anti-correlation ob- 
served by |Planck Collaboration! P011c[ ) (grey area shows the 
uncertainty). The black circles show the relation obtained after 
fitting the resulting SEDs in the Planck HFI and IRIS 100 jim 
bands for the central pixels of the series of filaments described 
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in Tab. 3, when the intrinsic (3- Tdust relation described by Eq 
is introduced in the Qats of LamC and aSil grains. In this case 
the modelled ficoiour - Tcoiour anti-correlation is also represented 
for the population of aggregates with a porosity degree of 40% 
(magenta squares) and the same intrinsic relation. 



dence is considered as a correction factor applied to a standard 
non-varying Q^bs 



Qabs(a,y) = eo(«,v)(v/y,)^«^""">^('"/^\ 



(8) 



where go is the grain absorption efficiency without the tem- 
perature dependence for size a and frequency v, v, - c/A, is 
the threshold frequency from which the correction applies, and 
SPiTdusi) = /3Q,,,{Tdus,) - Po is the index coiTection with re- 
spect to ySo, which is the spectral index without temperature 
dependence (Tab. 2). The threshold function in wavelength is 
H(x) - 0.5(1 + 4tanh(x)/5) where x - \og{AIA,) and s controls 
the steepness of the transition. We choose s - \ and /I, = 100 jum 
and deal with the clouds described in Tab. 3 for the DHGL pop- 
ulations. From laboratory measurements, [Mennella et al. |( |1998) l 
and Coupeau d et al.| ( |201 l| l demonstrated that for both amor- 
phous carbonaceous grains and amorphous silicates the opacity 
decreases with decreasing temperature, and that the spectral in- 
dex at a given wavelength increases when the grain temperature 
decreases from 300 K to 10 K. Hence we assume that for both 
LamC and aSil the emissivity spectral index varies with the tem- 
perature. For the sake of simplicity, we assume that the varia- 
tions are similar for the two dust populations with Po - 1 .55 for 
the LamC population and 2.11 for the aSil population. We do 
not consider the effects of noise (see Section [J!2] for a detailed 
description). 

Most of the available observations are more or less consis- 
tent with a Pcoiour - Tcoiour anti -Correlation that can be fitted 
with Pcoiour °^ ^OT/oi/; ' ^hcn considering the large error bars 
these fitted parameters suffer (fDesert et al. 2008 , Paradis et al. 
2010; Veneziani et al. 2010; Pla nck CoUabo ration 2 01 ld|c| ). For 
instance, Planck Collaboration! i 2011c| l observed starless cold 
clumps and found that 



^P 



Pcoiour = [(0.02 + 0.001) + (0.035 + 0.001)r,ow]" 

^-0.96 

-t ct 



'colour 



1.38 



(' colour \ 
20 K / 



(9) 
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for 7 < T colour S 20 K, when fitting the SEDs in the Planck HFI 
(857, 545, and 353 GHz) and IRAS 100 jum bands. To reproduce 
this relation, we have to introduce the intrinsic fi - T^ust relation 
in the absorption efficiency Qabs 



-0.96 



(10) 



which allows us to mimic 
15 < Toolour < 18 K (Fig. 



the observed anti-correlation for 
12 1. To get lower colour temper- 
atures, we apply the same intrinsic l3Q„,,XTdust) relation to the 
population of aggregates with 40% of porosity (j6o - 1.27) 
yielding colour temperatures as low as 12 K and almost the 
same ficoiour - TcoUmr relation as for the DHGL populations (see 
Tab. 2 and Fig. [T2] i. We repeated this exercis e in the Herschel 
PACS and SPIRE bandPjas in Para dis etal.N 2010), and in the 
BOOMERANG band^ as in |Veneziani et al.| ( [2010) . In the two 
cases, we were able to reproduce the observed anti-correlation 
using /3e^„,(r</„,,) ~ PcoiouriT colon,), independent of the spectral 
bands considered. 

The slope of the anti -correlation seems to be only marginally 
affected by radiative transfer effects for relations around ySg^,^^ oc 
Tt' This means that the observed relations reflect the intrinsic 

dust 

variations in the grain optical properties for SEDs fitted for the 
spectral bands of the instruments mentioned here when consid- 
ering noiseless observations. This might no longer be true when 
the relation is shallower, ySg,,,,, oc T^^^ -^ '° " as the radiative trans- 
fer flattening may then win over the intrinsic anti-correlation. We 
also find that the steepness of the relation is almost unchanged 
when changing dust populations from DHGL grains to aggre- 
gates and that the change in the spectral index without temper- 
ature dependence, ySo, from one dust population to another does 
not produce a strong discontinuity in the observed Pcoiour-T colour 
relation (Fig. 12 1. This indicates that the Pcoiour - Tcoiom rela- 
tions should be the same for all colour temperatures, at least as 
long as the clouds are starless (Mali nen et al.|[20TT| ), with the 
highest temperatures coming from the diffuse regions and the 
lowest from the centre of dense molecular clouds where grain 
growth is expected. We do believe that our simple modelling of 
the intrinsic ySg,,,,, - r^„s/ anti-correlation will be a very promis- 
ing ingredient of future studies including the detailed physics of 
the emission of amorphous solids (Meny et al. 2007), however 
this is beyond the scope of this paper. Finally, we note that the 
Pcoiour - T colour relations observed with Planck and Herschel are 
affected by noise, which may produce part of the anti-correlation 
(see Section[l2l Shetty et al. 2009, Blain et al. 2003). Therefore, 
the measured relations may overestimate the steepness of the 
true Pcoiour - T colour relation, where "true" means the relation ob- 
served after the radiative transfer effects but without the noise, 
and thus the steepness of the intrinsic Pg^iJJdusi) relation. 



6. Column density and dust emissivity estimates 

In addition to measurements of the dust colour temperature and 
emissivity spectral index, the dust SEDs from dense clouds can 
also be used to determine the hydrogen column density, Nh, or 



Paradis et al. 



j2010| actic plane at / = 59° with 

Her schel PACSMid SPIRE. 

Veneziani et al.'(2010) observed cirrus at high Galactic latitude 



(b ~ -40°) with BOOMERANG channels at 245 and 345 GHz 
(Balloon Observations Of Millimetric Extragalactic Radiation and 
Geophysics), and in the IRAS 100 /jm and DIRBE 240 /jm channels 
(Diffuse InfraRed Background Experiment onboard COBE, COsmic 
Background Explorer). 



the dust emissivity at a given wavelength, e{A). Knowledge of 
the column density is important as it allows us to estimate the 
mass of a cloud if its distance is known. The variations in the 
dust emissivity from one region to another is then commonly 
used as a tracer of grain growth from the diffuse medium to the 
centre of dense molecular clouds. In this section, we examine the 
accuracy of the estimates of Nh and e{X) using the dust SEDs, 
and also the influence of the spectral bands used to carry out this 
kind of analysis. 

6.1. Column density 

Using the colour temperatures determined with the fits of the 
dust SEDs, the column densities of the clouds can be expressed 



N„ = 



Sy 



By(Tcolour)KvfimH 



(11) 



where yu = 1.33 is the mean atomic mass, mn is the mass of 
the hydrogen atom, Sy is the emerging intensity at frequency 
V, By{T) is the Planck function at temperature T and frequency 
V, and Ky = e/jU»JH is the dust opacity at frequency v. We 
calculate the observed column densities for SEDs fitted in the 
Herschel PACS and SPIRE bands, and then in the Planck HFI 
and IRIS 100 yum bands (see Section [3T| . For the Planck esti- 
mate of Nh, we use the band at 857 GHz, and for the Herschel 
counterpart the band at 250 jum. The results for the ten clouds 
described in Tab. 3, containing DHGL grain populations, are 
shown in Fig. 13 In this case, the opacity at 250 /vm is equal 
to 0.051 cmVg and at 857 GHz to 0.025 cm-/g (Tab. 2). The 
column density is strongly underestimated with both Planck and 
Herschel bands, and this discrepancy increases with increasing 
Av = 1 to 20 from a factor of 1.8 to 3.0 in the Planck case, 
and from 2.5 to 6.6 in the Herschel case. This is in good agree- 
ment with previous studies (Evans et al. 2001; Stamatellos & 
|Whitworfli|2003) |Fischera|201 1, .Malinen etal..2011) . This dif- 
ference comes from the assumption of a single temperature in 
Eq. 1 1 The colour temperature is always higher than the tem- 



perature of the grains at the centre of the clouds and this over- 
estimate increases with the cloud central Ay as can be seen in 
Fig. [8] The influence of the difference between the colour tem- 
perature and the temperature of the grains on the calculation of 
Nh comes from the Planck function: for instance, at 857 GHz 
we found that B(18 K)/B(14 K) ~ 2, and at 250 /im that 
B(18 K)/B(14 K) ~ 2.6, which explains why the estimate of 
Nh with the Planck 857 GHz band is slightly more accurate than 
with the Herschel 250 /im band. These ratios increase to 9.5 and 
22, respectively, when the temperature differs from 8 K to 14 K. 
For comparison with more diffuse media and to test the influ- 
ence of the radiative transfer, we modelled a cloud with a central 
density equal to 30 H/cm-*, leading to a central visual extinction 
Ay - 0.09. The column density is underestimated by a factor of 
1.7 for Planck and 2.2 for Herschel. In this case, the discrepancy 
is not explained by radiative transfer effects but only by the mix- 
ture of several dust populations distributed in size and having 
different temperatures. 

Linear fits of the Nh - Ay relations corresponding to the dif- 
ferent populations of grains used in this paper, and for the cloud 
geometry and density distributions described in Se ction|2.3| are 
listed in Tab. 5. This is given for comparison with BohUn et al 
(1978) {Nh = 1.87 X lO^Uy H/cm^). Finally, real estimates of 
the column density from far IR and submillimetre observations 
are usually worsened by the ignorance of the true dust opacity 
inside dense molecular clouds. 
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Fig. 13. Case of the DHGL populations, a) True hydrogen col- 
umn density Nh at the centre of the modelled clouds as a func- 
tion of the central visual extinction Ay (red line). The dashed line 
shows the Nh -Ay relation measured by Bohlin et al. ( 1978| l for 
the diffuse ISM. The blue line represents the column density es- 
timated from the SED measured in the Planck-HFI and IRIS 100 
jum bands, whereas the black line comes from the SED-fits made 
in the Herschel bands, b) Ratio of the true column density of the 
modelled clouds to the column densities measured in Herschel 
and Planck bands (black and blue, respectively). 

Table 5. Relation between the hydrogen column density and the 
visual extinction for the dust populations Usted in Tab. 2. 
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6.2. Dust emissivity 

As we know the true column densities of the modelled clouds, 
we can measure the dust emissivity using Eq. |7] As in the pre- 
vious section, using the ten clouds in Tab. 3, we perform this 
analysis with the SED fitting in both Planck and Herschel bands. 



Figure 14 presents the variations in the dust emissivity at 250 
jL/m as a function of the visual extinction for the DHGL pop- 
ulations and for aggregates with a porosity degree of 40%. In 
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Fig. 14. Case of the DHGL grains (solid lines) and of the aggre- 
gates with 25% of porosity (dashed lines). Dust emissivity at 250 
fim as a function of the central visual extinction Ay. The black 
and the blue lines show the emissivity measured in the Herschel 
and in the Planck-HFI bands respectively for the central pixels 
of the filaments described in Tab. 3, assuming the observer does 
know the true column densities of the clouds. 



creases when the cloud density increases. This is due to the over- 
estimate of the dust temperature made when using the colour 
temperature, which was already seen in Section 6.1 for the mea- 



an opposite way to the variations observed in the ISM (Juvela 
|et aL||201l"| [Planck Collaboration! |2011d| l, the emissivity de- 



surement of Nh. The use of the grain temperatures at the centre 
of the clouds, as discussed in Section [43] would lead to an al- 
most flat relation. This shows that the increase in the submillime- 
tre emissivity towards dense clouds can be explained by neither 
purely radiative transfer effects nor the bias introduced by the 
use of the colour temperature. It must be caused by variations in 
the grain properties from the diffuse to the dense medium. Grain 
growth may explain the higher emissivity towards dense regions, 
as shown in Fig. 14 in terms of the variations in the emissivity 
of aggregates with 25% of porosity. For Ay = 10, the emissivity 
is increased by a factor of 4.8 compared with that of the DHGL 
grains. This factor is comparable to that observed in dense clouds 
with similar central Ay (see Tab. 1, Juvela et al. 201 1, Stepnik et 
al. 2003). 



7. Conclusion 

The primary aims of this paper were to differentiate between 
the effects of the radiative transfer and dust properties in dense 
molecular clouds. Secondly we have explored some explanations 
of the dust emission features observed towards these regions. To 
do so, we modelled the dust emission from idealized filaments, 
which were described as infinite circular cylinders at equilib- 
rium that were externally heated by the isotropic ISRF, includ- 
ing full radiative transfer calculations. We also briefly consid- 
ered the case of Bonnor-Ebert spheres (Appendix [B). Our main 
results are as follows: 

1. Radiative transfer effects, i.e. the mixing of different tem- 
peratures along the line-of-sight, can explain neither the low 
colour temperatures measured at the centre of dense and star- 
less molecular clouds, nor the observed ficoiour - T colour anti- 
correlation. The effect of radiative transfer is to produce an 
artificial correlation (not an anti-correlation as observed) be- 
tween the emissivity spectral index, Pcoiour, and the colour 
temperature. 
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2. For externally heated clouds with simple geometries (cylin- 
ders and spheres) containing dust populations with standard 
optical properties, the ficoiour - T colour relations observed with 
Planck and Herschel are too steep to be explained solely by 
noise. For starless interstellar clouds, it is thus likely that the 
observed anti-correlation is related to intrinsic variations in 
the grain optical properties with temperature as suggested 
by laboratory experiments ( Mennella et al. 1998 , Coupeaud, 
let al. 201 1| and mo dels of the emission of amorphous solids 
^Meny et al.|2007 1. However, an anti-correlation might also 
be produced if the sample includes map pixels with strong 
internal sources or (even a small number of) clouds observed 
with very low S/N's. 

3. The increase in the submillimetre emissivity from the dif- 
fuse to the dense medium cannot be explained by radiative 
transfer effects, but must originate in variations in the grain 
optical properties, which could be caused by their growth 
inside dense molecular clouds to form porous aggregates 
( |Ossenkopf & Henning|1994| ). 

4. Regarding grain growth, only the aggregation of grains in 
mixed aggregates of small and large grains can explain the 
deficit in emission in the mid-lR, the low colour tempera- 
ture, and the higher submillimetre emissivity observed to- 
wards dense molecular clouds. 

5. We found that for clouds with a central visual extinction less 
than 10 magnitudes, or central densities less than 2 x 10^ 
H/cm^ in the case of DHGL populations, a linear relation be- 
tween the colour temperature and the "true" dust temperature 
at the centre of the clouds does exist. Assuming that a reli- 
able measure of the cloud central visual extinction/density is 
available, it is possible to estimate the dust temperature using 
this relation. 

6. We investigated the impact of intrinsic variations in the grain 
optical properties with wavelength, I3{A). Even if possible, 
we showed that it should be rather difficult to track back the 
nature of the dust grains from the spectral variations in the 
dust SEDs emerging from dense clouds. This difficulty is due 
mostly to the radiative transfer effects for A < 300 fim and 
the mixture of different grain populations at longer wave- 
lengths. In real observations, the situation is further confused 
by the presence of noise. 

7. We also investigated the impact of intrinsic variations in the 
grain optical properties with temperature, /3(T). For the three 
observational results considered in this paper, with ficoiour ~ 
'^colour' '•^^ observed anti-correlation seems to reflect the 
/?- T^ust intrinsic relation. This result may no longer hold for 
shallower relations for which the radiative transfer flattening 
effect would win over the intrinsic anti-correlation. As noise 
may produce an anti-correlation, the observed relations may 
overestimate the steepness of the true jicoiour - Tcoiom rela- 
tions. 

8. Finally, we confirmed that the column density is strongly un- 
derestimated when determined by SED fitting. In diffuse me- 
dia, it is reduced by a factor of 1.7 and 2.2 in Planck and 
Herschel bands, respectively. This is explained by the mix- 
ture of several dust populations widely distributed in both 
size and temperature. This discrepancy then increases when 
the central visual extinction of the clouds increases because 
of the diflFerence between the colour temperature and the 
"true" dust temperature in the central layers of the clouds, 
where most of the material lies. 
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Appendix A: Details regarding the spectral 
variations in the emissivity spectral index 

We use the ten clouds described in Tab. 3 (diffuse medium case), 
illuminated by the ISRF, and consider five different mixtures of 
grains: LamC and aSil (see Tab. 2), ACAR and aSil, BE and aSil, 
LamC and ENSTA, and LamC and OLIVA (see Section|5?T]for a 
description of their optical properties). The corresponding SEDs 
and (Pcoimir - 7"co/oMr)-parameters, obtained with fits between 100 



and 3 000 yt/m, are shown in Figs. A.l i) and A. 2 



All the mixtures containing LamC grains have similar colour 
temperatures, Tcoiour, regardless of the silicate population (aSil, 
ENSTA, or OLIVA). If the LamC grains are replaced by ACAR 
or BE grains, the colour temperature drops. This decrease is 
due to their higher emissivity in the far-IR and shows that the 
colour temperature is dominated by the temperature of the car- 
bonaceous populations, which are always the warmest. 

The populations of ENSTA and OLIVA grains have very 
different optical properties in the submillimetre (Fig. This 
reappears in the fitted values of the emissivity spectral index, 
Pcoiour, which are much higher for a mixture containing OLIVA 
grains (1.75 < PcoUmr ^ 2) than for a mixture containing ENSTA 
grains (1.15 ^ (3 colour < 1-35, see Fig.[A!2i. 

Finally, Figs.|A.l 



) and c) show the spectral variations in the 
emissivities of the different mixtures of grains, which are mea- 
sured using Eq.|7] For all the mixtures containing aSil grains, for 
A > 500 /im, the spectral variations in the emissivity perfectly 
reflect the spectral variations in the absorption efficiencies, Qabs^ 
of the carbonaceous grains. For the mixture containing LamC 
and OLIVA grains, the carbonaceous grains still dominate the 
spectral variations in the emissivity for /I > 1 000 ytzm. However, 
if we replace the OLIVA by ENSTA grains, this is no longer 
true and the population of silicates dominates the spectral varia- 
tions for /I ^ 1 000 jum. This illustrates the strong dependence of 
the emissivity spectral variations on the dust populations consid- 
ered. However, the relatively simple behaviours described here 
are all mixed up if we consider that the optical properties of both 
the silicates and the carbons undergo spectral variations in the 
far-IR and in the submillimetre (see Section [5T| for a detailed 
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ACAR + aSil 
BE + aSil 
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10-' 10" 
Wavelength (pm) 

Fig. A.l. a) SEDs for the central pixel of the cloud with cen- 
tral Ay - 10 for a mixture of LamC and aSil (black), LamC 
and ENSTA (magenta), LamC and OLIVA (blue), ACAR and 
aSil (red), BE and aSil (green), b) Emissivity normalised at 100 
//m for the cloud with central Ay = 10, for a mixture of LamC 
and aSil (black solid line), of ACAR and aSil (red solid line), 
and of BE and aSil (green solid line). Absorption efficiencies 
normalised at 100 yum: aSil (black dotted line), LamC (black 
dashed line), ACAR (red dashed line), and BE (green dashed 
line), c) Emissivity normalised at 100 fim for the cloud with 
central Ay - 10, for a mixture of LamC and aSil (black solid 
line), of LamC and ENSTA (magenta solid line), and of LamC 
and OLIVA (blue solid line). Absorption efficiencies normalised 
at 100 /im; aSil (black dotted line), LamC (black dashed line), 
ENSTA (magenta dashed line), and OLIVA (blue dashed hne). 



description). As a result, it might be quite difficult to track back 
the nature of the grains from the spectral variations in the dust 
SED. 
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Fig. A.2. Emissivity spectral index as a function of colour tem- 
perature for fits of the dust SEDs between 100 and 3 000 
for the central pixels of the modelled clouds. Linestyles are the 
same as in Fig.lAk). 



Appendix B: The case of Bonnor-Ebert spheres 

To test the influence of the shape of the clouds on our results, we 
study the case of externally heated Bonnor-Ebert spheres, which 
are isothermal spheres at hydrostatic equilibrium (Bonnor 1956] 
|Ebert|[T955] l. We consider almost critically stable spheres with 
( = 6.5, or similarly pdpiR) - 14.1, where ^ is the dimen- 
sionless length equal to {AnGpc/c^y^^r and cs is the isother- 
mal sound speed, G the gravitational constant, and r the radius. 
We assume a constant kinetic gas temperature of Tg„s = 12 K. 

has a minor 



Fischera (201 1^ showed that the choice of f and Tgas • 
mtluence on the dust temperature, which depends mainly on the 
column density. The modelled clouds are divided into 50 con- 



centric cells (see Section 2.2 1 and the emission maps produced 
have 85 X 85 pixels with a single pixel size equal to 0.6% of the 
radius of the cloud. We vary the masses of the spheres between 
0.2 and 100 Mq to get central hydrogen column densities rang- 
ing from 2.5 x 10^" to 5 x 10^^ H/cm"^, or central Ay - 0.1 
to 17.7 in the case of DHGL grain populations. The parameters 
describing the spheres are shown in Tab. B.l. 

For comparison with the cylindrical clouds presented in the 
main body of the paper, we measure the ficoiour and Tciour pa- 
rameters for Bonnor-Ebert spheres with DHGL dust populations 
illuminated by the ISRF (Fig. B. 1 1, or by the ISRF extinguished 



by an external radiation field A"' up to 5 (Fig. 



B.2i. Figure B.3 
ummated by the 



shows the case where Bonnor-Ebert spheres, il 
ISRF, contain accC and accSi populations or aggregates with 
various porosity fractions. The values of the colour temperatures 
and the emissivity spectral indices vary when compared to the 
case of cylindrical clouds, although in all cases, the trends are 
similar. This means that the results described in the main body 
of the text in the case of cylinders hold for externally heated 
Bonnor-Ebert spheres. 

We now study the effects of noise on the Pcoiom — Tcoiom re- 
lations measured in Bonnor-Ebert spheres for observations in 
Herschel channels (PACS and SPIRE). The procedure is exactly 



the same as in Section 3.2 but we examine three cases depend- 
ing on the 250 yum peak surface brightness after convolution 
with the beam (the column density range covered by the mod- 
elled Bonnor-Ebert spheres is indeed larger than for the cylin- 
ders). For case A, we exclude the clouds with S/N at 250 yum 
lower than 8.3 (o S'aso/jm < 10 MJy/sr, excludes all clouds with 



M > 20 Mq), for case B lower than 4.2 (o 5250/im < 5 MJy/sr, 
M > 50 Mq), and for case C lower than 1.7 (o 5250/im < 2 



MJy/sr, M > 100 Mq). The results are presented in Fig. |B.4[ In 
case A, the correlation between Pcoiow and TcoUmr remains pos- 
itive. In case B, a slight anti-correlation appears {B > -0.6), 
but cannot account alone for the steep anti-correlation observed 
towards interstellar clouds {B < -1). In case C, where we in- 
clude the nine modelled clouds with the lowest S/N's, an anti- 
correlation appears and is steeper than B - -\ for more than 
half of the simulated cloud samples. The few faint clouds with a 
low S/N are responsible for the strong anti -correlation measured. 
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Fig. B.l. a) Same as Fig.[l|; in the case of Bonnor-Ebert spheres, 
b) Same as Fig.[TJl in the case of Bonnor-Ebert spheres. 
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Table B.l. Parameters describing the gas distribution of the Bonnor-Ebert spheres. M is the the cloud mass in units of solar masses, 
Nh is the hydrogen column density at the centre of the cloud, pc is the central hydrogen density, R is the external radius of the 
sphere, and Ay is the visual extinction at the centre of the cloud, which depends on the grain populations considered (labels are the 
same as in Tab. 2). 
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Fig. B.2. Same as Fig. [2] in the case of Bonnor-Ebert spheres. 
The blue line shows a sphere with Ay = 1.2 at the centre, green 
shows 6.3, red 12.5, and magenta 20.9. 
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Fig.B.3. Same as Fig.]?^ in the case of Bonnor-Ebert spheres. 



Fig.B.4. Same as Fig. |4j; in the case of Bonnor-Ebert spheres, 
for a signal-to-noise ration at 250 fj.m greater than a) 8.3 (case 
A), b) 4.2 (case B), and c) 1.7 (case C). 
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